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ABSTRACT 

We use the new Type la supernovae (SNe la) discovered by the Sloan Digital Sky Survey-II Su- 
pernova Survey together with additional supernova datasets as well as observations of the cosmic 
microwave background and baryon acoustic oscillations to constrain cosmological models. This com- 
plements the standard cosmology analysis presented by Kessler et al. (2009) in that we discuss and 
rank a number of the most popular non-standard cosmology scenarios. When this combined dataset is 
analyzed using the MLCS2k2 light-curve fitter, we find that more exotic models for cosmic acceleration 
provide a better fit to the data than the ACDM model. For example, the flat Dvali-Gabadadze-Porrati 
model is ranked higher by our information-criteria tests than the standard model with a flat universe 
and a cosmological constant. When the supernova dataset is instead analyzed using the SALT-II 
light-curve fltter, the standard cosmological-constant model fares best. This investigation of how sen- 
sitive cosmological model selection is to assumptions about, and within, the light-curve fitters thereby 
highlights the need for an improved understanding of these unresolved systematic effects. Our inves- 
tigation also includes inhomogeneous Lemaitre-Tolman-Bondi (LTB) models. While our LTB models 
can be made to fit the supernova data as well as any other model, the extra parameters they require 
are not supported by our information-criteria analysis. Finally, wc explore more model-independent 
ways to investigate the cosmic expansion based on this new dataset. 

Subject headings: cosmology: observations — supernovae : general 
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1. INTRODUCTION 

The Type la supernova (SN la) measurements which 
first indicat ed an accelerating expansion of t he universe 
(iRiess et al..,f 998; .Perlmutter et ah 19991 : see iFilippenkol 
l2005l for a review) have been confirmed and substan- 
tiated by a secon d generation of high-redshift super- 
nova experiments (iRiess et al.l 12004 : lAstier et al.l [20061 : 
IWood-Vasev et al.l 120071 : IRiess et al 1 120071). An impor- 
tant step forward in this respect was recently taken with 
the Slo an Digital Sky Sur vey-II (SDSS-II) Supernova 
Survey (jFrieman et al.|[2"008f ). This three- year survey, un- 
dertaken with a large CC D camera on a ded icated 2.5-m 
telescope in New Mexico (|Gunn et al.ll2006[ ). has discov- 
ered and followed several hundred SNe la, mainly in the 
redshift interval z = [0.01, 0.45]. These intermediate 
redshifts were previously underexplored, and filling this 
"redshift desert" not only provides impor tant new con- 
straints on cosmology (jKessler et al.|[2009t ). but will also 
help constrain systematic effects by bridging the low- 2; 
and the high-z supernova populations. 

The first-year SDSS supernova dataset is dis- 
cussed in three com panion papers, including this one. 
iKessler et al.l (|2009l hereafter K09) present the dataset 
in detail an d also use it to constra in standard cosmologi- 
cal models. iLampeitl et al.l ()2009f ) combine the SDSS SN 
data with other constraints to derive joint constraints on 
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dark energy from low-redshift (z < 0.4) measurements 
only; they also explore the consistency of the SN and 
BAO distance scales. 

To complement and extend the cosmological analysis 
presented in these papers we will in this paper use the 
first-year SDSS-II supernova data to explore several al- 
ternative c osmological models . Following the analysis 
outUned bv lDavis et al] (|2007l . hereafter DOT), we com- 
bine the 103 new SDSS-II SNe la in the K09 dataset 
with new analyses of previously available SN la datasets, 
as well as complementary data, to explore non-standard 
cosmologies. We also investigate the use of more model- 
free approaches in constraining the evolution of the uni- 
verse. 

This paper is organized as follows. In Section 2 we 
describe the datasets invoked in the analysis and how 
they are combined, while in Section 3 we present the 
cosmological models that are investigated in this work. 
In Section 4 we discuss our results. Section 5 includes a 
discussion of systematic effects, while Section 6 presents 
some ways of expressing generalized parameters from the 
supernova dataset. Finally, in Section 7 we provide a 
summary of our results. 

When we refer to the "standard model" we mean 
Friedmann-Robertson- Walker cosmology with a constant 
dark energy equation-of-state parameter, also known 
as "w Cold Dark Matter" (wCDM), of which the 
cosmological-constant model (ACDM) is a special case. 
We use units in which c = 1 . 

2. DATASETS 

In this paper we make use of constraints from several 
different datasets in order to test a number of cosmo- 
logical models. Compared to the previous analysis by 
D07, here we make use of an enlarged and re-analyzed 
supernova set, a new prescription as well as new data 
for the baryonic acoustic oscillations (BAO), and also an 
updated cosmic microwave background (CMB) analysis. 
In this section we present these datasets, and describe 
how they are combined. 

2.1. Type la Supernovae 

The primary new dataset used in this analysis com- 
prises the 103 ne w SNe la from the first-year SDSS-II 
supernova survey (jFrieman et al.ll2008t ISako et al.ll2008l : 
K09). 

This sample is published and discussed at length in 
our companion paper (K09) which also includes a com- 
prehensive and consistent re-analysis of other sets of lo- 
cal and high-z SNe la, using the same light-curve fitter. 
This is important since it ensures that all the supernova 
datasets are treated in a uniform manner regarding selec- 
tion criteria and light-curve fitting. For the analysis pre- 
sented in this paper we start by discussing the supernova 
dataset analyzed using the Multicolor Light Curve Shape 
2k2 fitter: (MFCS; Jha et al. 2007)^^. According to the 
MFCS analysis of K09, the standard ACDM cosmological 
model provides a rather poor fit to these data. We will 
also discuss calculations for the K09 S N dataset analyze d 
with the SALT-II light-curve fitter ijGuv et all I2007D . 

Wc have used the entire 

Nearby+SDSS+ESSENCE-I-SNLS-I-HST dataset (e) o f K09; 
|http://das. sdss.org/va/SNcosmology/sncosm09_fits. tar. gz| 



which is better fit by the ACDM concordance model. In 
total, we use 288 SNe from the analysis of K09. The dis- 
tance moduli and redshifts for these SNe are provided in 
K09 (their Table 10 for MFCS and Table 14 for SAFT-II). 
The very detailed and restrictive selection criteria for the 
total sample of SNe are described in K09; following their 
analysis we have used only supernovae with z > 0.02 and 
added an additional "intrinsic" dispersion of (Tadd = 0.16 
mag to the uncertainties output by the light-curve fit- 
ter. Since this added dispersion is motivated in K09 
to make the local MFCS SN Hubble diagram have a 
equal to unity per degree of freedom, it is essentially in- 
dependent of cosmology. A similar intrinsic dispersion 
(we use 0.14 mag) gives & Xdof—^ global SAFT- 

II fit. Our supernova dataset is t hus a well-selected 
sample compiled from rn a ny survevs (iKessler et al.ll2009l: 
Wood-Vasev et al.l 120071 : lAstier et al.l l2006l : iRiess et al.l 



2007t iJha et~ .1 120071 . and references therein) 



2.2. Cosmic Microwave Background 

When analyzing CMB observations there are two use- 
ful parameters commonly employed. One describes the 
scaled distance to recombination, TZ, and the other the 
angul ar scale of the sound horizon at recombination, £a 
(e.g.. iKomatsu et all 120091 lElgarOv fc Multamakil l2007t 
iWang fc Mukheriedl2006i r 

The shift parameter, TZ, is given by 




Ho dz 



H{z) 



(1) 



where Sk{x) — sinx', x, sinhx for k — 1,0,-1, respec- 
tively, and z* is the redshift of the last-scattering surface. 

The position of the first CMB power-spectrum peak, 
which represents the angular scale of the sound horizon 
at recombination, is given by 



_dA(f*) 

fs(^*) 



(2) 



where dp^{z^) is the comoving angular-diameter distance 
to recombination while the comoving sound horizon at 
photon decoupling, Tg, is given by 



H{z) 



-dz. 



(3) 



which depends upon the speed of sound before recombi- 
nation, Cs- Using both these parameters in combination 
reproduces closely the fi t from th e full CMB power spec- 
trum (but see also Elga rOv fc Mu ltamaki 2007, for some 
caveats), and within the standard model the two param- 
eters are only weakly correlated. Here we use the recent 
CMB measurements from the five-year Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) observations, and 
adopt the values I a = 302.10±0.86 and 11 = 1.710±0.019 
with c orrelation coefficient 0.1109 from IKomatsu et al.l 
(f2009f l. We further assume = 1090 exactly (Komatsu 
et al. 2009; variations within the uncertainties about this 
value do not give significant differences in the results). 

In a previous paper (D07), we used only the CMQ-TZ 
parameter to constrain cosmologies. The same method is 
used in the K09 paper. Here we will, for reasons outlined 
below, adopt a different approach. 
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Although TZ has been commonly used to constrain non- 
standard models, this approach may not always be en- 
tirely appropriate, because parameters close to standard 
loCDM were assu med in deriving the value of TZ (see, 
e.g.. Section 5.4.1 of lKomatsu et al]l2009l : [Kowalski et al.l 
l2008l Section 6). We will therefore use H. only for the 
wCDM model for which it was derived. The resulting fit 
is of interest because it can reveal any tension between 
the BAO/CMB and the supernova constraints. As will 
be further discussed below, ACDM is not a good fit to 
the data when MLCS is used for the SN analysis, and this 
motivates our search for better fits among non-standard 
models. 

We therefore instead perform an analysis without TZ for 
all the models. In doing this, wc use only the product of 
the acoustic scale £a with the position of the BAO peak 
(Section 12. 3p to complement the SN data. 



0.0 



2.3. Baryon Acoustic Oscillations 

As with the CMB, there are several parameters in com- 
mon use for comparing BAO observations to theoretical 
models. The most immediately observable of these is a 
measurement of the ratio of the sound horizon scale at 
the drag epoch, rs{zd), to the dilation scale, Dv{z). The 
drag epoch, Zd ~ 1020, is the epoch at which the acoustic 
oscillations are frozen in. 

A more model-independent constraint can be achieved 
by multiplying the BAO measurement of rs{zd)/ Dv{z) 
with the CMB measurement £a = '^dA{z*)/rs{z:t:), thus 
cancelling some of the dependence on the sound horizon 
s cale. 

iPercival et al.l (|2009l ) measured rs{zd)/Dv{z) at two 
redshifts, z = 0.2 and z = 0.35, finding rs{zd) / Dv{0.2) = 
0.1905 ±0.0061 and r^^Zd) / Dv{0.ib) = 0.1097 ± 0.0036. 
Combining this with £a gives the combined BAO/CMB- 
£a constraints: 



dA{z*) rs{zd) 



DviO.2) rsiz,) 
dA{z*) rs{zd) 



DviO.35) rsiz,) 



18.32 ±0.59, 



10.55 ±0.35. 



(4) 



This combin ation is equiva l ent to the Sk/Dy combi- 
nation used in iPercival et al.l (j2007f ). but with the ratio 
of the sound horizon at the two epochs made explicit. 
Before matching to cosmological models we also need to 
implement the correction for the difference between the 
sound horizon at the end of the drag epoch, Zd « 1020, 
and the sound horizon at last-scattering, z* w 1090. The 
first is relevant for the BAO, the second for the CMB, 
and rs(zd)/rs(Zi,) = 1.044 ± 0.019 (using values from 
iKomatsu et al.l [20091 ). Inserting this into Equation |4| 
gives the final constraints we use for the cosmology anal- 
ysis: 



dAjz*) 
DviO.2) 
dA{z*) 



17.55 ±0.65, 



10.10 ±0.38. 



(5) 



i:»v(0.35) 

We take into account the correlation between these mea- 
surements using the corr e lation coefficient of 0.337 cal- 
culated bv IPercival et all ()2009[ ). 
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Fig. 1. — Flat dark-energy model (Fw): a flat universe with 
constant w. The constraint from each of the observational probes is 
shown by shaded contours. These are all 95% confidence intervals 
for two parameters. Overlaid with black lines (95% and 99.9% 
confidence intervals) are contours from combining CMB/BAO-^a, 
CMB-7^ and SN constraints. The shaded contour labeled SN is 
for the analysis using the MLCS light-curve fitter. In this plot we 
have also added the CMB-7?. constraints, although these are not 
included in the model selection. The dotted supernova contours 
arc using the SALT-II fits. For the SALT-II dataset the combined 
contours are given by the dashed contours, and are clearly in better 
agreement with the cosmological-constant value, w = —I, shown 
by the dashed-dotted line. 

The ratio of sound horizon distances between drag 
epoch and last-scattering was calculated using the FA 
model (see Table [1]) for the evolution between those two 
redshifts. Wc expect this to be a good approximation for 
all the models we test here because the redshift difference 
between the decoupling and the drag epoch is relatively 
small, and the sound horizon at decoupling and drag is 
mostly governed by the fractional difference between the 
number of photons and baryons. 

2.4. Combining the Datasets 

To clarify how we combine the data, we show in Fig- 
ure[l]our best fit to the data in the w — ^lm plane for a flat 
universe with constant w (The Fw model, see Table |T]) . 
The constraints from each of the observational probes arc 
shown by contours (according to the figure legend). In 
this and in all similar figures, these are 95% confidence 
intervals for two parameters. The supernova dataset fit- 
ted with MLCS is shown by the shaded (red) contours, 
while the SALT-II SN analysis is shown by the dotted 
contours. The combined contours (95% and 99.9% con- 
fidence intervals) are overlaid in black for MLCS and 
by the dashed contours for SALT-II. In our calculations, 
we are marginalizing over a common magnitude shift for 
all SNe, thus allowing for arbitrary values of the Hub- 
ble constant, Hq, and the typical absolute magnitude of 
SNe la. 

This standard cosmology case is the only one for which 
we also show the CMB-7?. constraint. The constraints 
from the prescription for the BAO that we have imple- 
mented, which also include the CMB-^^i parameter, are 
labeled CMB/BAO in the figures. Using this CMB/BAO 
product cancels out some of the dependence on the sound 
horizon size at last scattering. This thereby removes the 
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dependence on much of the complex pre-recombination 
physics that is needed to determine that horizon scale. 

In Figure [1] the likelihood analysis takes into account 
the correlation between £a and TZ. In other words, our 
likelihood plot shows contours derived from CMB-7<^ sep- 
arate fr-om contours derived from BAO/CMB-£^ that are 
then combined taking the weak correlation into account. 
For the other models, and for the model selection, the 
CMB-7?. is not included. It is then the combination of 
BAO / CMB-£a with SN data that we use as the basis for 
the model comparison performed in this paper. This new 
analysis is more model independent, and therefore better 
suited for ranking exotic models. Excluding the TZ con- 
straint also leaves more room for models with non-zero 
curvature. 

It was noted by iPercival et al.l ()2007f ) that the dis- 
tance between the two BAO redshift bins seemed to be 
in slight conflict with the supernova measurements for 
the standard model. Compared to the Percival et al. 
(2007) BAO/CMB constraints this tension was exacer- 
bated by our new MLCS SN dataset. The new BAO 
data ease some of this ten sion (see also IPercival et al.l 
120091 : iLampeitl et al.l l2009f ) . Rather than selecting only 
one of the BAO redshift measurements, or any one of the 
SN light-curve fitters for that matter, we choose to ex- 
plore the implications of such tension for the more exotic 
models of interest in this paper. 

Figure [T] also helps us understand in a more quanti- 
tative way how the different datasets constrain the cos- 
mology. The combined fit using MLCS for the SNe (the 
SALT-II case will be discussed further in Section 14. 3p 
and also CMB-7e gives w = -0.79 ±0.13, while the 
more general method we employ for all models in this 
paper (omitting TZ) gives a best fit at w = —0.83 ± 0.24. 
The error bars quoted are 95% confidence level for one 
parameter, and we only consider statistical errors. For 
the approach used by K09 (SN+CMB-7e+BA0-A) we 
find w = —0.78 ±0.13. Although the procedure we 
have chosen to combine the different observational con- 
straints also affects the results, it is clear that it is the 
MLCS analysis of the supernova dataset that makes the 
main difference in pushing w away from the cosmological- 
constant value of w = —1. We therefore move on to 
explore more exotic cosmological alternatives. 

3. TESTING NON-STANDARD MODELS 

iKessler et all (I2009D introduced the first-year SDSS-II 
SNe to supplement existing data in order to constrain 
the standard cosmological model. They concentrated on 
testing a cosmological-constant model and a flat universe 
model with constant equation-of-state parameter. Here 
we will extend that analysis to a more general investiga- 
tion including several of the most popular non-standard 
models. The models are briefly presented below. 

3.1. Beyond- Einstein Models 

There are many specific models based on new funda- 
mental physics that make predictions for the expansion 
history of the universe. We follow D07 in using infor- 
mation criteria to rank these models; see als o the r ecent 
investigations by, e.g., iKurek fc Szvdlowskil (j2008f ) and 
iRubin et al.l ([20091). 

With the new SDSS-II supernova data and updated 
data on the first-peak position of the CMB, combined 



TABLE 1 



Summary of models 


Model 


Abbrev. 


Parameters^ 


Flat cosmo. const. 


r A 




Cosmological const. 


A 


O 

1 £m , i 


Flat constant w 


Flo 




Constant w 


w 


f^m, fife, W 


Flat w(a) 


Fwa 


flm, WQ, Wa 


Cardassian 


Ca 


^m, q, n 


Flat Chaplygin 


FCh 


A 


Chaplygin 


Ch 


Qfc, A 


Flat Gen. Chaplygin 


FGCh 


A, a 


Gen. Chaplygin 


GCh 


Qky A, a 


DGP 


DGP 




Flat DGP 


FDGP 




LTB Gauss 


LTBg 


f2in, Qout, ro 


LTB Sharp 


LTBs 





The free parameters in each model. When 
fitting the SN la data we also fit an additional 
parameter, M, for the normalization of SN mag- 
nitudes. We include this in the number of degrees 
of freedom and in the number of free parameters 
considered when calculating information criteria, 
but do not list it here as a parameter in each 
model. For more details of these models and the 
parameters, see D07. 



with the new approach to the BAO constraints, we show 
below that it is worthwhile retes ting the non-stand ard 
cosmological models examined bv lDavis et al.l ()2007t ). 

The selected models include both exotic dark-energy 
models as well as alternative models that can be in- 
terpreted in terms of modifying the theory of gravity. 
The first class of models we test are the standard 
cosmological constant, constant w, and variable w 
models. We also include Dvali-Gabadadze-Porrati 
(DGP), Cardassian expansion, and several versions of 
Chaplygin gas models. The appropriate references and 
the equations we use for describing H { z) in e ach of these 
models were collected by iDavis et al.l ()2007L Equations 
7-18), and we refer the reader to that paper for further 
information. For reference, in Table [T] we list all the 
different models included in this study. 



3.2. The Lemaitre-Tolman-Bondi Models 

In addition to the above models, we also test some 
inhomogencous cosmologies. Such models have gath- 
ered significant interest in recent years as a means to 
explain the cosmological observations without invoking 
dark energy. In the simplest class of such models we 
live close to the center of a large, spherically sym- 
metric void d escribed by the Lemaitre-Tolman-Bondi 
(LTB ) metric ()Lemaitrelll93l Il997t Ilblmanlll934t iBondil 
Il947f ). The apparent acceleration of the expansion 
is then caused by the spatial gradients in the met- 
ric, such that our local region has a larger Hubble pa- 
rameter than the outer region. While the LTB mod- 
els challenge the Copernican principle, several stud- 
ies have shown that they cannot be ruled out by 



Enavist & MattssonI 12007 


: iGarcia-Bellido & HauEb0llel 


2008: Caldwcfl & Stebbins 


I2008D. 



The LTB models are characterized by two arbitrary 
functions, often expressed as the expansion rate H{r,t) 
and the matter density parameter Om(?', i), which de- 
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pcnd not only on time, but also on the radial coordinate. 
As a consequence, inhomogeneities arise independently 
in the matter distribution and in the expansion rate. 

We will consider LTB models constrained by the re- 
quirement that the Big Bang occurred simultaneously 
throughout space by implementing a particular choice 
of i?oW, 



Ho{r) 



1 



^ujr) 



sinh 



-1 



I Mr) 



(6) 



so that the time of the Big Bang was tsB ~ i^o 
for all observers irrespective of their position in space. 
The model is then completely specified by only one 
free function, the matter density parameter flMir), with 
^M{r) + rifc(r) = 1. We consider two different density 
profiles. In the first model, i^uir) takes the form of a 
Gaussian undcrdensity. 



(7) 



This model has three free parameters, where il-m is the 
matter density at the center of the void, fiout is the 
asymptotic value of the matter density, and ro is the 
scale size of the underdensity. 

In the second model, fluir) has a much sharper tran- 
sition from the local to the asymptotic value. 



1 + e"''"/'^'' 

1 + g{r-ro)/Ar 



(8) 



Here tq characterizes the size at which the transition oc- 
curs and the extra parameter Ar characterizes the tran- 
sition width. In the limit where Ar goes to zero, the 
density profile becomes a step function. The two models 
are designated LTBg and LTBs in Table [TJ 

Our sh arp LT Bs model is equivalent to the constrained 
model in lGarcia-BcUido fc Ha ugboUcj (|2008l ). We fix the 
transition width to be Ar = 0.065 Gpc h~ to obtain a 
sharp transition, so this is not a free parameter in our 
model. Furthermore, in both our LTB models we allow 
r^iii and Slout to take any positive values < 1, i.e., the 
models need not be asymptotically flat and in principle 
we allow also for solutions with a local overdensity. 

4. RESULTS 
4.1. Model Ranking Using the MLCS Analysis 

The results of all our fits using the MLCS-fit SNc and 
the combined CMB/BAO constraints are summarized in 
Tabic [2 The corresponding results for the SALT-II-fit 
SNc arc presented in Section 14.31 The results arc stated 
in terms of and the given degrees of freedom, as good- 
ness of fit (GoF), and in terms of information-criteria 
(IC) assessments. 

The backg round for the use of IC for these models was 
reviewed by iDavis et all (|2007l, their Sec tion 2). We 
follow their approach (see also lLiddl3l2004f ) and use two 
IC to select the best-fit models. These are the Bayesian 
information criterion (BIG) 

BIG = -21n/: + fclniV, (9) 

and the Akaikc information criterion (AIC) 

AIG = -21n/: + 2A:, (10) 
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Fig. 2.— The constraints from SNc and CMB/BAO on the 
parameters in the DGP model. The results have changed substan- 
tially from those of Davis et al. (2007). This is both due to the new 
datasets and our choice of not using CMB-7^. The flat DGP model 
is indicated by the vertical dashed-dotted line; for the MLCS fit, 
it is the best-ranked model by the IC analysis. The SALT-II fit 
to the SNe is again shown by the dotted contours. The combined 
constraints using the SALT-II SNc outlined by the dashed contours 
represent a poorer match to the CMB/BAO for the flat model. 



where C is the maximum likelihood, k is the number of 
parameters, and N is the number of data points used 
in the fit. A ABIC larger than 6 would be considered 
unsupported as compared to the best model. 

The number of degrees of freedom is derived from the 
288 SNe and the two CMB/BAO measurements, less 
the contribution from M and the number of parame- 
ters listed in Table [TJ Note that the fairly low per 
degree of freedom stems from the fact that we add an in- 
trinsic dispersion to the supernova data. As mentioned 
above, this is derived by measuring the dispersion in the 
nearby sample. The SDSS SNe actually have a lower 
dispersion (0.08 mag), which K09 attribute to selection 
effects (K09, their Appendix E). We also ran a test with 
this lower intrinsic dispersion, 0.08 mag, and found that 
while the fits arc now much worse (higher and very 
low GoF), the relative A AIC are not affected by much. 
Some rank order differences would be seen in Tabic 2, 
but only for differences of AIC<1, and these should not 
be regarded as significant. Since the intrinsic dispersion 
has somewhat different meaning in the SALT and the 
MLCS frameworks, it makes little sense to compare the 
GoF for the different light-curve fitters. 

Compared to the analysis performed earlier by D07, 
we can see that the new MLCS dataset provides different 
results in several interesting ways, as follows. 

• There are now a number of exotic models that 
fare rather well under the IC test. It is clear that 
both this dataset and the method we have used 
to combine the observational constraints permit 
a larger variety of models. We note in particular 
that many of the models provide very similar x^ 
values, and the IC tests therefore are sensitive to 
the number of parameters for the given model. 



6 



SoUerman, MortscU, Davis, Blomqvist, et al. 




0.4 0.5 0.6 0.7 0.8 0.9 

A 



Fig. 3. — Constraints for the flat generalized Chaplygin gas 
(FGCh). From the constraints in this analysis, this model cannot 
be excluded. The dashed-dottod a = line corresponds to param- 
eters that match the A model (with Qm = 1 — ^), and these do not 
match the combination with the MLCS fit data. The best MLCS 
fit is at o = -0.4, A = 0.55, which is far from the a = 0.0, A = 0.7 
that represents the best fit to the FA model. This is different from 
the analysis of DOT when the best fits were acquired for parameter 
values that mimic the cosmological constant. The combined con- 
straints from the SALT-II SNe (dashed) are more consistent with 
FA values. 

• The simple, fiat, cosmological-constant model (FA) 
favored by D07 is no longer on top of the ranking 
list in Table [H This can of course be understood 
from Figure [T] Instead, the model favored by both 
ABIC and AAIC is the flat DGP model that was 
unsupported by previous studies. In Figure [2] we 
can see the constraints on the DGP model from the 
new data and analysis. These are very different 
from the case presented by D07. The change is 
driven by the new SN data, and is exacerbated by 
the way we now combine the CMB/BAO and SN 
constraints (i.e. omitting CMQ-TZ). 

• The flat dark-energy model (Fif) may still be a vi- 
able model for this new dataset in terms of AIC. 
However, in Figure [T] we see that the most likely 
value from the combined dataset is w = — 0.83 ± 
0.24 (95% confidence level [C.L.] for one parame- 
ter). For this model, the value of fim ~ 0.31 ±0.06. 
Including CMB-7?, as in Figure [1] which is viable 
for this model, gives w ~ —0.79 ± 0.13. This off- 
set from ui = — 1 is clearly a feature of the new 
MLCS analysis of the combined supernova dataset 
(see also K09). The more general Ywa model gives 
the best-fit value of wq = —0.73 ± 0.38, while there 
are no useful constraints on the time varying com- 
ponent of 'w{a). 

• Models such as the generalized Chaplygin gas were 
found by D07 to be good fits to the data only 
when their parameters mimicked a cosmological 
constant. The new best fits instead fall in regions 
of parameter space that do not correspond to the 
cosmological constant. This is, for example, shown 
for the flat generalized Chaplygin gas in Figure [3l 

The Cardassian expansion model is not very well 



TABLE 2 

Information-Criteria Results for MLCS 



Model 


X I dof 


GoF (%) 


AAIC 


ABIC 


AAAICiii, 


FDGP 


2SS ri / 288 


99.21 


0.0 


0.0 


—2.3 


\ 


23S 2 / 287 


99.12 


1.9 


5.6 


0.0 


DGP 


233.2 / 287 


99.12 


2.0 


5.6 


0.4 


FGCh 


233.3 / 287 


99.11 


2.1 


5.7 


0.3 


Ch 


233.4 / 287 


99.10 


2.1 


5.8 


-1.7 


Fw 


233.5 / 287 


99.09 


2.2 


5.9 


0.5 


Fwa 


233.1 / 286 


99.03 


3.8 


11.1 


1.9 


w 


233.2 / 286 


99.02 


3.9 


11.2 


1.0 


Ca 


233.2 / 286 


99.02 


3.9 


11.2 


0.7 


GCh 


233.2 / 286 


99.01 


3.9 


11.3 


0.9 


FA 


237.3 / 288 


98.70 


4.0 


4.0 


4.0 


LTBg 


235.5 / 286 


98.69 


6.2 


13.6 


4.2 


LTBs 


237.6 / 286 


98.31 


8.3 


15.7 


6.8 


FCh 


257.6 / 288 


90.09 


24.3 


24.3 


13.2 



Note. — For our SN sample analyzed with MLCS plus 
CMB/BAO, the FDGP model is preferred by both the AIC and 
the BIG. The AAIC and ABIC values for all other models in the 
table are then measured with respect to these lowest values. The 
goodness of fit (GoF) approximates the probability of finding a 
worse fit to the data. The models are given in order of increasing 
AAIC. The final AAAIC^i, column simply displays the difference 
in ranking when using different prior s in the MLCS supernova 
analysis, as described in Section l4.1.1l 

constrained but also in this case the best fit does 
not match the FA model. This agrees with the 
poorer rating of the cosmological-constant model 
in the present investigation. 

• The only model that is very strongly unsupported 
compared to the other models is the flat Chaply- 
gin gas. 

In addition to the contour plots, we provide in Figured] 
a Hubble diagram displaying a selection of the discussed 
models. This representation offers perhaps a less abstract 
way to judge the models. Note that the best-fit models 
are in the upper panel of Figure[3]constrained both by the 
SNe la and the CMB/BAO, i.e., the model parameters 
are not optimized for a supernova Hubble diagram. 

4.1.1. Changing the MLCS Pnors 

The MLCS fight-curve fitter used by K09 is an updated 
version of the MLCS used bv IWood-Vasev et all (|2007l . 
hereafter WV07) . The new features include in particular 
new priors and an updated treatment of extinction. As 
shown by K09, this new treatment significantly changes 
the best-fit w value for the Fw model, and this is not a 
feature of the new SDSS-II SN dataset. 

To illustrate how these parameters within MLCS affect 
our cosmology ranking we have also performed calcula- 
tions for an alternative set of assumptions. Following the 
discussion in K09 (their Section 10.1.4 and Fig. 32) we 
start with the WV07 prior and then implement all of the 
cumulative changes up to (but not including) the point 
where K09 also change Ry from the typical Milky Way 
value of 3.1 to the lower value of Ry = 2.2. This is of 
course a rather arbitrary choice of parameter changes, 
and is just meant to illustrate how different MLCS as- 
sumptions and implementations affect the fits to the 
Hubble diagram for the models tested in this paper. This 
particular choice corresponded to an offset in w of < 0.1 
for the Nearby+ESSENCE+SNLS dataset analyzed us- 
ing the Fw model according to K09 (their Fig. 32). We 
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Redshlft 

Fig. 4. — (Upper) Hubble diagram for the MLCS supernova 
analysis. Distance modulus differences in magnitudes are shown 
with respect to an empty universe. The grey points are the SN la 
distance moduli with error bars, while the black points arc binned 
data. The binning is done using nA2=constant, where n is the 
number of SNe in the bin and Az the red shlft range. This uniform 
and unbiased binning was introduced by IRiess et al.l ((2007) to en- 
able visualizing the Hubble diagram with bins distributed over the 
entire redshlft range, and with nAz = 6 we get 8 bins with error 
bars that are of similar size in both directions. The redshlft error 
bars show the standard deviation of the redshifts in the bin, while 
the distance modulus error bars give the standard deviation of the 
mean of the distance modulus uncertainties within the bin. Note 
that the binning is only used for visualization in this figure, and 
that all the calculations and fits are performed on the unbinncd 
data. We have also plotted the best fits for a sample of models 
(we chose mainly the non-flat models in order not to clutter the 
figure). In the upper panel, these best fits are using both MLCS 
SNe and the CMB/BAO constraints, i.e., they are not optimized 
for the SN Hubble diagram alone. In order to plot these Hubble 
diagrams, the best fit value has been used for the different mod- 
els. For example, for the upper panel and the A model we derived 
M-M = Slog (c/Hq) + 25 = 43.35. Here, c is the speed of light, 
Hq the Hubble constant (where (c/Hq) is in Mpc) and M is the 
absolute magnitude of a Type la supernova. The other models 
included in the figure have similar values for M (c = 0.02 mag), 
with the LTBs model having the largest offset. Note that in all the 
calculations, these nuisance parameters were marginalized over. 
(Lower) Hubble diagram for the best fits to the MLCS SN la data 
alone, relative to an empty universe. This panel is zoomed in to fo- 
cus on the sharp LTB (dotted) and the Gaussian LTB (dot-dashed) 
models. For comparison we have also plotted the best-fit A model. 



perform here the calculations for all of our exotic models, 
for our choice of CMB/BAO prior, and for the complete 
SN dataset used in our analysis. 

The results of the different assumptions in the MLCS 
are given as the AAAIC/ju values in the last column of 
Table [5] This is simply the differences in AAIC com- 
pared to our original MLCS analysis, i.e., AAAICi^u — 
AAICmlcs ~ A AlCnew priors, with both runs normalized 
to the best-fit model. To not overload the presenta- 
tion, we have chosen to include only the AIC statistics 
here, since this is an asympto tically unbiased estimator 
(jBurnhamn fc AndersonI f2004f ) . The differences are not 
very large in this representation. Note again that dif- 
ferences of <1 in AIC are not significant. All models 
fare slightly better compared to the FDGP model in this 
case. A noticeable difference is that the FA model fares 
significantly better using this set of MLCS priors. In 
the Fw model, w = —0.99 ± 0.28. We stress that we in 
no way prefer this set of parameters, but use them only 
to illustrate how such changes would affect our model 
selection. 

4.2. LTB 

Here we take the opportunity to further discuss the 
results of our LTB models. Combining the SN data 
with the CMB/BAO constraint provides a critical test 
of the validity of our LTB models, since there appears 
to be some tension between the best-fit parameter val- 
ues preferred using the SNe and CMB/BAO data sep- 
arately. This of course assumes that the BAO con- 
straints can be used in this way also fo r LTB mod- 
els (see iGarcia-BeUido fc Haugb0lIi 120081) . The total 
values for our LTB fits are comparable to those of 
the other models (Table [2]), and because of the ex- 
tra parameters the LTB models fare poorly in the IC 
tests. These conclusions genera lly agree with what 
iGarcia-Bellido fc HaugboIIg (|2008f ) found for their LTB 
models. 

Our best-fit models (using the MLCS supernova anal- 
ysis) show the transition, rg, at 1.2 (1.8) Gpc for the 
LTBg (LTBs) models. For the LTBg the matter density 
f^iii = 0.25, r^out = 0.72, and these numbers are similar 
for the LTBs model. This is illustrated in Figure [5] 

We note that some of the LTB models studied in 
the literature display distinct features in the low- z 
end of the Hubble diagram (e.g., Alexan der et al.ll2007l : 
iChfton et al.][200l . It was therefore anticipated that the 
SDSS-II data would provide decisive constraints by fill- 
ing in the redshifts desert between z « 0.1 and z « 0.4. 
The smooth Hubble diagram in this redshlft regime pro- 
vided by our well-calibrated SDSS-II supernova dataset 
(see K09, their Fig. 23) clearly disfavors such hmrinosity 
distance vs. redshlft relations. 

However, not all LTB models investigate d pos- 
sess such con spicuous features (jAlnes etall [20061 : 
iGarcia-Bellidofc Haugb0llc 200j). In fact, our simple 
LTB models can fit the SN data (alone) with a some- 
what lower than that of the FA model. This is illustrated 
in the Hubble diagrams in FigureHl where the LTB mod- 
els are included together with a selection of other alterna- 
tive models. In the upper panel Hubble diagram all mod- 
els have also been constrained by the CMB/BAO data. 
This combination produces a poorer fit for the LTB mod- 
els and pushes the transition to higher redshifts. The 
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Fig. 5. — The LTB Gauss model is able to give a reasonable fit 
to the data but is not favored by the IC tests. It is clear that the 
data prefer the matter density to vary from a low value locally to 
a higher value in the distant universe. The scale size of the void in 
the best-fit model is 1.2 Gpc/h.. The SALT-II SNe (dotted) push 
f2in to lower values while leaving the constraint on flout unchanged. 

lower panel instead displays the LTB Hubble diagram 
when fit only to the supernova data, and demonstrates 
that a good fit to the data can be obtained. 

The more general conclusion of our LTB model test is 
therefore that the SDSS-II SN set is useful for testing 
the LTB models, but perhaps not decisive to the degree 
previously anticipated. Although the SDSS-II supernova 
set strongly constrains LTB models that display egre- 
gious features in the intermediate-redshift Hubble dia- 
gram, versions of the LTB model that do not display 
such features, such as the LTBg and LTBs models con- 
sidered here, can be made to fit the supernova data. 

However, in terms of our IC test, we have found that 
these LTB models are not well motivated; the number of 
extra parameters involved is not justified by a better fit 
to the SN data, in particular not when the CMB/BAO 
constraints are also included. Note that these results also 
hold when using the SALT-II light-curve fitter to analyze 
the SN data (Table [3]) . We also mention that if we were 
to impose the prior of asymptotic fiatness on our LTB 
models, the for these fits would increase even further. 

4.3. Model Ranking Using the SALT-II Analysis 

We have also performed calculations using SN data fit 
wi th the comp letely independent light-curve fitter SALT- 
II (|Guv et al.|[2007il , following the implementation as de- 
scribed by K09. K09 noted a substantial offset for the 
value of w as calculated in the standard model using 
MLCS versus SALT-II, and devote a long discussion to 
the differences of these light-curve fitters. As seen in Fig- 
ure [l] the SALT-II fits are more consistent with the con- 
cordance model, and we here investigate in detail how 
such a difference alters our ranking of the investigated 
exotic models. 

In performing these calculations, we had to exclude 
four of the SNe that were well fit by MLCS because in 
the SALT-II analysis they gave poor fits^^. 

24 The SNe designated d085, e020, k429, and HAWK. 



The SALT-II prescription clearly makes the SN data 
more compatible with the standard models. In Table [3] 
we see that the highest-ranked model is the FA model, 
a fiat universe with a cosmological constant. It has a 
similar as many of the other models, and is primarily 
favored by the IC since it has only one free parameter. 
This echoes the conclusions of D07. In the rightmost 
column of Table [3] we provide the difference in AAIC 
compared to the MLCS results, to illustrate which mod- 
els benefit most from using the SALT-II SN set. Note 
that the numbers of degrees of freedom are less than for 
the MLCS table both since we had to exclude four SNe, 
and due to the two extra free parameters in the SALT-II 
fits. We note the following: 

• FA evidently fares much better with SALT-II than 
with MLCS, as was already evident from Fig- 
ure [TJ For the Fw model the best-fit value for the 
equation-of-state parameter is w = —1.14 ± 0.26 
(95% C.L. for one parameter), with fl^ — 0.30 ± 
0.05. Including CMB-7?, as in Figure [1] which is 
viable for this model, gives w ~ —0.98 ± 0.14. 

• On the contrary, the FDGP model is now pushed 
toward the bottom of the ranking list. Figure [2] 
again illustrates the situation, with the SALT- 
II SN constraint dragging the solution along the 
CMB/BAO constraint and away from the fiat uni- 
verse. The non-fiat DGP model is still one of the 
models that fares better in this investigation than 
in D07 (as is the Ch model). 

• The LTB models fare as poorly (or worse) com- 
pared to the best-ranked models when the SALT-II 
SN data are used. In this respect the conclusions 
from Section 14.21 hold independent of the adopted 
light-curve fitter. 

• The SN distances obtained using the SALT-II light- 
curve fitter provide results that are overall more 
consistent with the analysis of D07. For the FGCh 
model (Figure [3]), we can see that the SALT-II SN 
distances allow a match to the CMB/BAO con- 
straints for parameter values that are more con- 
sistent with those expected for the cosmological- 
constant model. 

5. SYSTEMATICS 

When performing the same kind of IC ranking of cos- 
mological models as in D07, the new datasets and meth- 
ods explored in this paper lead us to quite different con- 
clusions. When we use the MLCS analysis the standard 
cosmological-constant model no longer stands out as the 
best and simplest model, and new room is made for a 
range of more exotic models. The results from the SALT- 
II analysis are quite different. 

A basic underlying assumption in the analysis is that 
the errors are mainly statistical. The systematic errors 
in, for example, t he SN survey s are e x tensively discussed 
elsewhere (e.g., lAstier et al.l 120061 : IWood-Vasev et all 
l2007t iKessler et ani2009[ ). and the uncertainties include, 
in particular, the distribution of host-galaxy extinctions 
(Ay) and the dust-reddening properties (Ry)- 
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While model selection techniques are not well equipped 
to face systematic errors, we can regard the above- 
mentioned calculations as a simplistic perturbation the- 
ory approach to test how assumptions regarding light- 
curve fitters affect the model selection. 

In Table [2] we illustrate how the model selection is af- 
fected by rather arbitrarily changing the priors in the 
MLCS light-curve fitter. Given the change in results en- 
gendered by using the ratio of total-to-selective absorp- 
tion of Rv = 2.2 that K09 derived from the SN sample 
itself, such a low value of Ry should clearly find support 
also from independent astrophysical observations. 

The most dramatic differences noted in this investiga- 
tion come from the very different results obtained for the 
two light-curve fitters investigated. This clearly points to 
unresolved systematic differences that need to be further 
explored. Table [3] shows in the final column that the fits 
of some exotic models are dramatically affected by the 
choice of light-curve fitter in the a nalysis of the SN data . 

While we acknowledge the way iKowalski etaD (|2008D 
treat the SN data from their large UnionOS compilation, 
we also note that adding a dispersion to the error mim- 
ics the use of more free p arameters in the fi t. Using 
the same UnionOS dataset, iRubin et al.l (|2009f ) reached 
conclusions similar to the ones in this paper; the differ- 
ence in values for different cosmological models are 
rather small. We have not included any systematic er- 
rors in our tests, and including larger systematic errors 
would clearly make it even more difficult to differentiate 
between cosmological models for a given dataset. 

However, it is also clear from our investigation that 
trying to encapsulate the systematic errors is notori- 
ously difficult. The substantial disagreement in model 
selection obtained with different light-curve fitters in- 
vestigated here demonstrates that deeper understanding 
of the systematics is warranted. It is also important to 
refrain from simply choosing the method or datasets 
that confirm previous findings or present prejudices. 

6. GENERAL EXPANSION HISTORY 

Given the difficulty in assessing particular models, we 
also consider some ways to derive more general cosmolog- 
ical results from our dataset. The usual parameter-fitting 
techniques used in cosmology are limited by the neces- 
sity of assuming a model with parameters to fit; there is 
more information in the cosmological data than can be 
extracted by parameter fitting. In this section we will 
ignore specific dark energy theories and instead try to 
elucidate directly the expansion history of the universe, 
H{z), or the evolution of an equation-of-state parameter 
of dark energy, w{z). The SN dataset fitted with MLCS 
is used here, since this easily provides distance moduli 
without having to fit a specified cosmology for the H{z) 
section, and it avoids additional free parameters for the 
w{z) analysis. 

6.1. The Equation-of-State Parameter w{z) 

We start this investigation by assuming a piecewise 
constant equation-of-state parameter, in order to fit 
a more general dark-energy model. The value Wi is 
calculated in each redshift bin i (defined by the up- 
per limit Zi), where we have arbitrarily chosen Zi = 
[0.15, 0.6, 1.7, 1090], roughly corresponding to low-, mid-, 



TABLE 3 

Information-Criteria Results for SALT-II 



Model 


XV dof 


GoF (%) 


AAIC 


ABIC 


AAAICmlcs 


FA 


274.4 / 282 


61.63 


0.0 


0.0 


4.0 


FGCh 


273.0 / 281 


62.32 


0.6 


4.3 


1.5 


A 


273.4 / 281 


61.66 


1.0 


4.6 


0.9 


DGP 


273.4 / 281 


61.65 


1.0 


4.7 


1.0 


Ch 


273.9 / 281 


60.87 


1.5 


5.1 


0.7 


Fw 


274.0 / 281 


60.67 


1.6 


5.3 


0.6 


GCh 


273.0 / 280 


60.61 


2.6 


10.0 


1.3 


Fwa 


273.4 / 280 


60.05 


3.0 


10.3 


0.8 


Ca 


273.4 / 280 


59.97 


3.0 


10.3 


0.9 


w 


273.5 / 280 


59.82 


3.1 


10.4 


0.8 


LTBg 


276.5 / 280 


54.71 


6.2 


13.5 


0.0 


FCh 


281.9 / 282 


49.06 


7.5 


7.5 


16.8 


FDGP 


282.5 / 282 


48.05 


8.1 


8.1 


-8.1 


LTBs 


289.3 / 280 


33.90 


18.9 


26.2 


-10.6 


Note. 


— For our SN sample analyzed with 


SALT-II plus 



CMB/BAO, the FA model is preferred by both the AIC and the 
BIG. The AAIG and ABIC values for all other models in the table are 
then measured with respect to these lowest values. The goodness of fit 
(GoF) approximates the probability of finding a worse fit to the data. 
The models are given in order of increasing AAIC. The numbers of 
degrees of freedom are loss than for the MLCS fits, since four SNe were 
omitted and because SALT-II uses two extra fit parameters. The final 
column simply displays the AIG differences as compared to the MLCS 
results in Table[2] i.e., AAAICmlcs=AAICmlcs-AAICsalt- 

high-z SNe and CMB data, respectively. We have fitted 
both Wi and ilm using a Monte Carlo Markov Chain, 
assuming a fiat universe. 

Note that the Wi are correlated; the covariance matrix 
is not diagonal. We c an, however, decorrelate the Wi es- 
timates by (following iHuterer fc Cooravl l2005f) changing 
the basis through an orthogonal matrix rotation that di- 
agonalizes the covariance matrix. This corresponds to 
applying a weighting function to the Wi. The uncorre- 
lated Wi are thus linear combinations of all Wi described 
by the weight function. 

That is, appropriately weighted data from all redshifts 
are used to obtain w in each bin, and in our case, the 
value of Wi corresponding to low redshifts (wi ) is a linear 
combination of the value of w{z) at z < 0.15 (~ 82% 
contribution), w{0.15 < z < 0.6) (~ 18% contribution), 
and negligible contributions from higher redshifts. In the 
same way, W2 is a combination of w{z < 0.15) (~ 24% 
contribution), w(0.15 < z < 0.6) (~ 73% contribution) 
and w{0.6 < z < 1.7) 3% contribution). The value of 
W3 has a ~ 17% contribution from w(0.15 < z < 0.6) and 
a ~ 87% contribution from w(0.6 < z < 1.7), whereas 
W4 predominantly comes from the highest redshift bin 
and only has small negative contributions from the three 
lower redshift bins. 

The results are presented in Figure[6l and are discussed 
in Section [Ql 

6.2. The Hubble Parameter H{z) 

IWang fc Tegmard (|2005f ) proposed a technique to ex- 
tract the expansion history in uncorrelated redshift bins 
from S Ne la data. This technique was used to great ef- 
fect bv lRiess et aT1(|2007( ). We now apply it with the new 
information from the SDSS-II supernova sample. 

The results are displayed in Figure [71 where the theo- 
retical models have been normalized as described in the 
caption. The redshift bins have been chosen by the re- 
quirement that the ratio of the error on the Hubble pa- 
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Fig. 6. — Observational constraints on a pieccwisc constant to (z) 
plotted (diamonds) at the middle value in each redshift bin. Note 
here that results are given not in terms of the original w{z) but 
weighted with data from all redshift bins in order to decorrelate 
the bins. The weights are displayed as thin lines with styles cor- 
responding to the respective points, and with the scale given on 
the right-hand side of the figure. This is further detailed in the 
text (Section [5TTJ- The results include our CMB/BAO constraints. 
The error bars are 68%. The constraints on the evolution of w{z) 
at 2 > 1 are clearly very weak. 

rameter and the redshift range probed by each measure- 
ment should be close to a constant times the redshift 
derivative of the Hubble parameter; we want to keep a 
uniform signal-to- noise ratio across the bins. Ideally, this 
constant should be equal to one, corresponding to 



dH ^ 

dz 



(11) 



where an is the uncertainty in the Hubble parameter and 
is the redshi ft range that each measurement probes, 
as calculated by IWang fc Tegmar^ {2003) • We compute 
dH/dz using a fiducial cosmological model, here taken to 
be a flat universe with rim = 0.3. 

We find that when using three redshift bins, the error 
from the width of the redshift bin dominates, whereas 
for four bins, the Hubble parameter uncertainty is some- 
what larger than the corresponding error from the red- 
shift range probed. Note that only the supernova data 
are used in obtaining H{z). 

6.3. The Deceleration Parameter q 

Another less theory-dependent investigation comes 
from exploring the deceleration parameter and its evolu- 
tion. Following iMortseU fc ClarksonI ()2009D . we consider 
a flat universe where 

q{a) = qa + qi{l - a) = qo + qi— — • (12) 

1 + z 

At zero redshift, g(a = 1) = go a-nd in the infinite past, 
q{a = 0) = qo + qi- This parametrization appears to be 
reasonably flexible in the sense that performing a least- 
squares fit to the dark-energy models employed in this 
paper always provides an acceptable fit. The comoving 
distance is given by 



dc{z) 



H{z) 



1 

Ho 



TT / exp 



[1 + q{u)]du 



(l+u) 



dv , 
(13) 



Fig. 7. — With the rich SN dataset we can also trace H{z) in a 
more model-independent way. Here the H{z) evolution is compared 
to some theoretical models. The models are normalized to agree 
at redshift zero, and we have assumed Ho = 72 km s^^Mpc^^. 
The data points are normalized to have the lowest-redshift bin on 
the w = —1 curve, and the four different points are independent 
of each other. Only supernova data were used for this plot, and a 
flat universe was assumed. 




qo 



Fig. 8. — Another less model-dependent investigation comes 
from e xploring the deceleration pa rameter and its evolution. Fol- 
lowing [Mortseir&||Clarkso^ II2009I) . we plot this representation to- 
gether with our new observational constraints. On the abscissa, the 
evidence for current acceleration is evident. The dashed-dotted line 
demarcates models with early deceleration (above the line) from 
model with early acceleration; that the best-fit lies well above this 
line demonstrates evidence for past deceleration. 

and the luminosity distance, which is the relevant quan- 
tity for SN la observations, is given by 



1 



V\^\Hodciz) 



(14) 



The angular diameter distance, relevant for the scale of 
BAO and CMB, is given by dA = dL/{l + z)'^. We can 
now use Figure [H] to study the expansion history of the 
universe within this simple flat model, as discussed be- 
low. 

6.4. General Results 
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Fig. 9. — The growth of structure is an independent way to 
constrain cosmological models. This is illustrated h ere f or standard 
gravity (SM) and for the DGP model (see Section [63J . It is clear 
that the present data do not yet select between these models, but 
any future deviations from 7 = 0.55 above a few percent would 
indicate deviations from ordinary gravity. Solid and dashed lines 
are 95% confidence levels for one parameter while the filled yellow 
(light) and orange (dark) regions show 68.3% and 95% confidence 
levels for two parameters respectively. 

The method of presenting the evolution of w{z) and 
H{z) given in Figs. [6] and [7] may have more general and 
long-lasting value than fits to specific models. 

The H{z) evolution shown in Figure [7] is displayed to- 
gether with a few cosmological models. We can see how 
the DGP model produces an acceptable fit to the evolu- 
tion of the Hubble parameter. The general trend, an in- 
creasing value for H(z) with increasing redshift, is clear. 

Any w{z) evolution would be an obvious signa- 
ture of exotic models. We do not, however, see 
any statis tically significant evo lution of w in Figure [6] 
(see also iKowalski et al.l |2008[ ). This representation 
could also be used to constrain generic "freezing" an d 
"thawing" models, follow ing, e.g.. I Sullivan et all (|2007f ): 
ICaldwell fc Linden ()2007D , and in this respect Figure [S] 
can be seen as an illustration of our ability to model- 
independently constrain the evolution of w using cur- 
rent data. We note that the result that w is consistent 
with —1 across the redshift range does not depend on 
the choice of binning, although the size of the error bars 
does. Using Eq. 11 we derive for our combined SNc la 
and CMB/BAO constraints (Figure^ go = -0.34±0.16 
and qi = 0.93 ± 0.25. The negative value of qo demon- 
strates that the present universe is accelerating. Note, 
however, that with the current (MLCS) dataset, this con- 
clusion does not follow for this parametrization for the 
SN dataset alone. Furthermore, the fact that the com- 
bined contour clearly hovers above the dashed-dotted line 
demonstrates that the universe was also decelerating in 
the past. The redshift where the universe transits from 
a decelerated expansion to an accelerated expansion is 
constrained to 0.41 < zt < 0.72 (95% C.L. for one pa- 
rameter). 

6.5. Growth Factor, 7 

There are a number of more exotic models that cannot 
be ruled out, even in principle, if the magnitude-redshift 



evolution we observe remains consistent with ACDM. 
This is because they can mimic standard dark energy 
perfectly in this respect. 

For many of these models, it would be useful to investi- 
gate additional and complementary constraints that can 
distinguish them from ACDM. This can be done using in- 
formation abo ut structure formation as expressed by the 
growth factor (|Lindeij|2005l : iLinder fc Cahnll2007f ). It is 
a measurable parameter of a model, in the same sense as 
the matter density Qm or the equation-of-state parame- 
ter w. All models in which general relativity holds have 
a growth parameter 7 « 0.55, but models that explain 
the acceleration by modifying gravity may have a signif- 
icantly different growth factor. Thus, this approach has 
the potential to separate models that have identical pre- 
dictions for expansion history but differ in the form of 
gravity. 

We show in Figure [9] the predicted growth factor for 
the DGP model and for the standard model (SM=Fu', 
dashed and solid contours, respectively), together with 
constraints from current data (filled contours). The 
DGP model has the same number of parameters as the 
cosmological-constant model, and can therefore provide 
a useful framework for testing modified gravity alterna- 
tives. The DGP model fared very well in the IC test in 
Table 2, but does give different predictions for the growth 
of structure. 

Figure [H] uses the new data from iGuzzo et all (l2008f) 
analyz ed following the procedure of ILinder fc CahnI 
(|2007f ). The filled contours are growth-factor data (/ c>c 
using / = 0.49±0.14 at z = 0.15 and / = 0.91±0.36 
at z = 0.77. It is clear that the present data do not 
yet select between these models, but any deviations from 
7 = 0.55 above a few percent would indicate departures 
from ordinary gravity. 

7. CONCLUSIONS 

We have followed the information-criteria model selec- 
tion outlined in D07 using the new supernova dataset 
from the SDSS-II (K09) and a new way to combine the 
updated BAO and CMB datasets. The ranking using 
the MLCS supernova analysis is quite different from the 
models selected by D07, with more room for exotic mod- 
els. Using instead the SALT-II light-curve fitter, we ob- 
tained fits more consistent with the cosmological con- 
stant. Some LTB models were also tested; they were 
able to give decent fits to the SN data, but were not 
favored by the full IC tests. 

We also discuss how more general results can be de- 
rived from this new dataset, and demonstrate that the 
expansion of the universe is indeed speeding up now, 
while it used to be slowing down. This conclusion is 
based on combining constraints from the supernova data 
with complementary data. 

The information-criteria tests used in this analysis are 
fairly simple, but serve the purpose to quantitatively 
rank models with rather similar values. The way 

the results differ between the different light-curve fitters 
suggests that more detailed statistical tests are not jus- 
tified. 

The data from the releases of several major supernova 
searches (SDSS-II, ESSENCE, SNLS, HST) have been 
analyzed in this paper. Most of these teams plan to 
release substantial updates of their supernova samples 
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soon. This will be important in improving the statistics, 
but our analysis also makes clear that systematic effects, 
in particular uncertainties in the treatment of extinction 
and light-curve fitting, will limit the predictive power in 
selecting a model based on these data. Rather than de- 
termining w with increased accuracy, sor netimes claimed 
to be measured to within ~ 6% (e.g., iKowalski et al.l 
l2008t iKomatsu et al.|[2009[ ). this work cautions on the 
importance of the systematic effects. 

We note that the recent investigation by iHicken et al.l 
(|2009[ ) also finds differences in the derived value for w us- 
ing different light-curve fitters and also for different cuts 
in the datasets, and provide several good suggestions on 
how to improve future supernova cosmology in this re- 
gard. The SDSS-11 datasct is indeed of sufficiently high 
quality to enable further studies of the relevant systcm- 
atics. 

We have explored how systematics limit our ability to 
select the best models for the universal expansion based 
on current data. A discussion on the benefits of MLCS 
versus SALT-II is beyond the scope of this paper. K09 
trace part of the discrepancies to the way the light-curve 
fitters treat extinction, how this affects very blue SNe and 
in particular to uncertainties in the bluest pass-bands. 
Better constraints on the systematics from the final data 
releases from SDSS-II, ESSENCE, and SNLS may wcU 
help in resolving some of the tension discussed in this 
paper. We may then find that our standard model holds. 
Or we may simply live in a more exotic universe. 
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